On the MHD load and the MHD metage 
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In analogy with the load and the metage in hydrodynamics, we define magnetohydrodynamic 
load and magnetohydrodynamic metage in the case of magnetofluids. They can be used to write the 
magnetic field in MHD in Clebsch's form. We show how these two concepts can be utilised to derive 
the magnetic analogue of the Ertel's theorem and also, how in the presence of non-trivial topology 
of the magnetic field in the magnetofluid one may associate the linking number of the magnetic 
field lines with the invariant MHD loads. The paper illustrates that the symmetry translation of 
the MHD metage in the corresponding label space generates the conservation of cross helicity. 

PACS numbers: 47.65.-d, 52.30.Cv 
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I. INTRODUCTION 



In a pioneering paper [lj, Lynden-Bell and Katz intro- 
duced the concept of load and metage in the hydrody- 
namical flows to show that the content of Kelvin's circu- 
lation theorem is contained in the conservation of a par- 
ticular function of load; if that very function is same for 
two flows then the flows are isocirculational. In the case 
of MHD flows, one can define similar analogous quanti- 
ties that, as we shall see in this paper, can prove to be 
quite significant both physically and mathematically. 

A magnetofluid (also termed hydromagnetic fluid or 
MHD fluid) of infinite conductivity but possibly com- 
pressible and viscous is quite an interesting class of flu- 
ids in which the flow velocity (v) and the magnetic field 
(B) interact actively; in general, B cannot be consid- 
ered to be a passively advected quantity. Anyway, owing 
to the Alfven's theorem, B remains frozen in the mag- 
netofluid; as we shall see, this property helps to design 
three active scalars that remain conserved on the fluid 
particles allowing themselves to be used as the particle 
labels which by definition, although vary continuously 
throughout the fluid, remain fixed on the fluid parti- 
cles and hence serve as possible coordinates of the label 
space in the Lagrangian description of the flow. With the 
comment that these scalars may become passive scalars 
in certain well-studied cases e.g., in kinematic dynamo 
problem wherein one basically studies the instability 
around B = state of MHD because there the magnetic 
field itself is treated to be passively advected, let us now 
look for the aforementioned advected scalars. 
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II. MHD LOAD AND MHD METAGE 

We take an isolated volume of magnetofluid of negli- 
gible resistivity. B will be penetrating it all throughout. 
Due to solenoidal nature, some magnetic field lines will 
be entirely inside the magnetofluid and some will pene- 
trate out of the boundaries and complete the loop out- 
side the magnetofluid. Lets assume that the later type 
of field lines begin at some the initial point on the sur- 
face which we term as "entry" to leave the magnetofluid 
at some corresponding final points on the boundary that 
we call "exit" . Having taken such a line, we construct a 
magnetic flux tube of small cross-section around the line 
which serves as a sort of axis to the tube. Obviously, for 
the magnetic field loops lying entirely within the mag- 
netofluid, the "entry" and "exit" collapse onto the same 
point and hence one should keep in mind that in the cal- 
culation that follows one would need to introduce a "cut" 
somewhere on the flux tube defined about such a loop. 
We also assume for the time being that the magnetic field 
lines are not knotted. 

Let AS(l) be the cross-section of the tube at a distance 
/ from the entry and at that point let the mass of the in- 
finitesimal disk of fluid of thickness dl be dm. So, mass 
of the tube (up to the first order in AS{1)) is 



Am 



dm 



entry 



p(l)AS(l)dl 



(1) 



entry 



Again, if Atfi is the strength of the flux tube then similarly 
one has 



A4> = B(l)AS(l) 



(2) 



Using the relation © in the equation |T]) and taking the 
limit as the cross section of the flux tube going to zero, 
one arrives at 
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^ _ dm 



entry 



B(l) 



d! 



(3) 
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where A may be called "magnetohydrodynamic load" in 
analogy with the "load" in hydrodynamics. Now, we give 
this scalar A the status of the scalar field by defining A 
at a point as the A of the field line passing through that 
point. As, A is constant along the magnetic field line, so 
it must satisfy the equation: 

B.VA = (4) 

i.e., A = A is a solution of the relation ([4}; we have 

5.VA = (5) 

and as A remains constant with the motion of the field 
line due to the frozen condition of the magnetic field, 
we can write the following mathematical form for the 
condition: 



DX 
~Dt 



= 



(0) 



Again, let A = A be another solution of the relation (|3|) 
such that 



dX 
dX 







(7) 



So, one has like the equation ([5]), a similar equation for 
A that we shall call "conjugate magnetohydrodynamic 
load" due to the reasons which will become obvious in due 
course. A naturally satisfies (by definition) the following 
equation: 

B.VA = (8) 

From the equations ||SJ| and ([8]) we can write 

B = C*(VA x VA) (9) 

because B is orthogonal to both VA and VA. Here, due 
to the solenoidal nature of the magnetic field, one may 
conclude taking divergence of the relation © that 



C = C(A,A) 



(10) 



Redefining, A as J* C(A, A)c?A we can do away with C to 
get ultimately from the relation 



B = VA x VA 



(11) 



which is the form of magnetic field associated with Cleb- 
sch parametrised magnetic vector potential. Thus, now 
B automatically satisfies the relation V .B = 0. A 
very interesting geometrical interpretation for the rela- 
tion (fTTj) is that the part of the field line inside the mag- 
netofluid is basically the curve formed due to the inter- 
section of the surface of constant A with the surface of 
constant A. Thus the specification of the families of the 
such surfaces is an alternative way of specifying the B in 
the magnetofluid; A and A though do not vary along the 



field lines, distinguish between the field lines. 

A third scalar is now needed to complete the set of 
coordinates of the label space. Let us refer the equation 
([3]) again and change the upper limit from "exit" to some 
floating point (which we shall call "float"), on the narrow 
tube about the field line, at a distance I from the "entry" . 
Using the concepts of conservations of mass and magnetic 
strength, define the following Lagrangian invariant scalar 
which we define as "magnetohydrodynamic metage" ji: 



dmf 



float 



entry 



B(l) 



dl 



(12) 



where the subscript "/" on m refers to the fact that one 
is interested only in the mass of the tube up to the length 
ending at the floating point. Due to the frozen nature of 
the magnetic field line, the floating point moves with the 
fluid and this makes ^ a candidate for labeling the fluid 
elements. One, thus, has 



~~Dt 



= 



Also, from the relation (fP2|) it follows that 



B.Wfi 



(13) 



(14) 



From the equations 
the Jacobian J is 



(fTTjl and [fl4]). one finds that 



J 



d{x,y,z) 



(VA x VA).V^ 



(15) 



Now this value of J confirms that A, A and fi are the inde- 
pendent variables in the Lagrangian description and we 
can think of a label space with coordinates (A,A,/z). The 
coordinate system (A,A,/z) may be treated as curvilinear 
coordinate system in the usual (x,y,z) space. One can 
choose VA, VA and V/x as basis vectors in this curvilin- 
ear coordinate system. As A is a label so it must satisfy 
the following equation: 



DA 
~Dt 







(16) 



So just like A, A is a Clebsch type co-moving scalar field. 
Also, from the relations ([6]), (fl"6|) and (jlip the usual in- 
duction equation written below must follow. 



OB -» -* 
dt K ' 

v being the velocity field of the MHD flow. 



(17) 



III. A AND MAGNETIC ERTEL'S THEOREM 

One can always define a vector potential A for B. Al- 
though A is always non-unique in the sense that gradient 
of any scalar field may always be added to it without 
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changing the B and hence one would notice that in the 
arguments that follow it hardly matters which A one is 
choosing. Anyway, w.r.t. the basis vectors of the coordi- 
nate system (X,X,jS), A may be written as 



A = A1VA1 + A 2 VX 2 + A 3 VA 3 



(18) 



For convenience, we have renamed (A,A,/j,) as (Ai, A2, A3). 
Lets define 



B„ = Vx x A 



(19) 



where, Va = (8/8X1,8/8X2,8/8X3). Seemingly, the re- 
lation ([TO]) suggests that B p is the magnetic field mea- 
sured in the curvilinear coordinate system (A,A,/x). We 
shall see (equation (f26"|) ) that each component of B p is a 
conserved quantity which may be interpreted as a poten- 
tial magnetic field B\ w.r.t some Lagrangian invariant 
scalar. 

Potential magnetic field B\ w.r.t any Lagrangian in- 
variant scalar, say A for example, may defined as 



B x = 



B 



.V A 



(20) 



Due to the frozen nature of the magnetic field which is 
the result of the following equation valid in the MHD of 
infinite conductivity: 



(21) 



As the infinitesimal displacement 5x between two mov- 
ing fluid particles follow equation similar to the equation 
(|2~Tj) . we have 




dSx 
~dt 



= (Sx.V)v 



(22) 



Due to the relation one can easily put the equations 
(f2"Tj) and (|2"2"|) to good use and can churn out in the pro- 
cess the following equation: 








(23) 



that may be simply interpreted as the conservation of B\ 
(i.e. potential magnetic w.r.t. MHD load) on fluid parti- 
cles; obviously this is the analogue of Ertel's theorem[3j 
in MHD. 

The vector potential A (see relation (02)), may also be 
written as 



A = A[i + A' 2 i + A' 3 k 



(24) 



The transformation rule for the components of A in two 
different basis vectors (refer to relations (|4"2")l and (|2"4"|) is: 



(25) 



This rule ([25)1 and the relation (fT5[) allows one to establish 
that 



(V A x A) j = B Xi 



(26) 



This is what we wanted to arrive at. Comparing the 
relations (fl9|) and (|26|) we can see that the magnetic field 
measured in the curvilinear coordinate system (A, A, /it) is 
a conserved quantity. 



IV. MHD LOADS AND LINKING NUMBER 

A useful function of A may be defined in a steady MHD 
flow. Steady, by definition, means the explicit indepen- 
dence of the physically measurable quantities from time; 
hence p and B do not depend on time in the case of 
steady MHD flow. The definition of A (equation ©) 
would then suggest that the magnetohydrodynamic load 
as well is not dependent on time. Lets define the function 



$(A) = J B.dS 



(27) 



that basically gives the strength of the flux through a 
large magnetic field tube surfaced by a locus of points 
having constant A; for convenience of discussion we 
henceforth shall mean such a tube when we shall be using 
the phrase "tube of load A" and shall denote the mass 
within it by M(X). From the definitions (J3]) and (|2T|) . we 
can very well write 



dM _ d<$> 
~dX~~dX 



(28) 



If A = Ao is the MHD load at the center line of the tube 
of load A, one may note that 



M (Ao) = 0; $(A ) = 



(29) 



Hence, from the equation (|28p one notices that $(A) and 
M(X) and MHD load can be related to each other by dint 
of the following relations: 

M(X) = [ X X'^^ldX'= [* A($'W (30) 
Ao dX Jo 

xi dM ^ dx' = r w ™l (3D 



$(A) 



A„ A' dX> 



n 



A(M') 



One of the interesting mathematical applications of defin- 
ing <&(A) is that through it the concept of MHD load A 
can go on to express magnetic helicity in its own terms 
as should be obvious from what follows. 

Consider the case that the magnetic field lines spiral 
around toroids which lie entirely within the magnetofluid. 
Such an arrangement would inherently introduce mag- 
netic helicity in the flow. This should be evident from 
the fact that one can build such a structure by "Dehn 
surgery", i.e. by cutting a closed magnetic flux tube at 
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a section, twisting the free ends through a relative angle 
that is integral multiple of 27r, and reconnecting; for sim- 
plicity's sake one may suppose that the resulting twist 
is uniformly distributed round the tube. Obviously, two 
fluxes of magnetic lines of force are associated with each 
toroid: one being the flux threading the meridional cross- 
section of the toroid and the other being the one which 
penetrates the inner boundary of the equatorial section 
of the toroid. Let the strengths of the fluxes be $ m and 
$ e respectively. In this non-trivial topology, one can de- 
fine two invariant MHD loads to characterise the toroid 



as 



dM 



dM 



(32) 



where M is the mass inside the torus that provides 
the magnetic field lines a compact support V on whose 
boundary dV, n.B = 0. 

Now, a conserved pseudo-scalar quantity viz., magnetic 
hclicity[4j, |5j, ly] is given by: 



H 



A.BdV 



(33) 



The integrand here measures how much A rotates about 
itself times its modulus; thereby making the name mag- 
netic helicity apt because it gauges the relative curling 
and braiding of the lines of A and B. Thus, H gives an 
intuitive measure as to what degree the field lines resem- 
ble helices. The most important property of H which 
is of our concern is that it takes non-vanishing values 
only if the topology of the integral curves of magnetic 
field is non-trivial. Obviously in the definition (|3"3"|) . A 
is not unique, for, a term V0, 8 being a scalar field can 
always be added to it keeping B unchanged. Let us pon- 
der over this non-uniqueness of the vector potential Q. 
For smooth discussion's sake, we assume for the time 
being that the volume is simply connected. Suppose 
Ai — > Ai + di6, then from the definition (|3"3"j) of mag- 
netic helicity we can find the change 5H in H to be: 



SH 



ve.Bd A x 



6B.hd 2 x 



(34) 



dv 



where n is the unit vector perpendicular to the infinitesi- 
mal surface element d 2 x and we have used the solenoidal 
nature of magnetic field and Gauss divergence theorem. 
The relation (|34[) amounts to saying that the magnetic 
helicity will be gauge invariant in case the surface dV 
bounding V is a magnetic surface i.e., B.n = on dV . 
This condition for gauge invariance is rather strong be- 
cause \i B.n 7^ on dV then one cannot seek refuge in 
Coulomb gauge for it is too loosely defined inside V with 
no information about the outside field whatsoever. More 
starkly, it means that different solenoidal vector poten- 
tials inside V can correspond to Coulomb potentials of 
fields which have different structures outside V. Now if 
we relax the condition that V is simply connected, then 



the line integrals of A about the 'holes' in the possibly 
multiply connected region have to be specified in order 
to have gauge-invariant magnetic helicity within dV on 
which B.h = 0. 

It is well-knownQ that if the flux strength in the 
tube before the 'Dehn surgery' was $ and the number 
of twists required to get the toroidal configuration men- 
tioned above be n {i.e., the linking number of any pair 
of the magnetic field lines is n) then 



H = n& 



(35) 



But since in the configuration under study, one may use 
relations (|3"2"|) to rewrite the definition (|3"3")l as: 



H = 



M 



*e(M') 



dM 



7 $m(M') 



dM' 



(36) 



It is of interest to note that using the relations (|32p , (|33|) 
and (|36p . one may express the linking number in terms 
of the invariant MHD loads as in the following relation: 



A I 



M' 



dM"dM' 



1 



* 2 Jm>=oJm»=o Xe(M")X m (M') 



(37) 



This evidently provides a physical and topological signif- 
icance for the MHD loads. 



V. fi AND CROSS HELICITY 

Like the topological conserved quantity magnetic he- 
licity, in the non- viscous MHD with infinite conductivity 
one has yet another conserved quantity viz. , cross helicity 
defined as follows: 



H r = 



\BdV 



(38) 



It measures the degree of mutual knottedness[5[ of the 
vorticity field and the magnetic field; this remains con- 
stant even though the vortex lines are no longer frozen 
in the fluid. In this section, we shall prove that an ap- 
propriate symmetry displacement associated with the in- 
finitesimal change in the MHD metage [i generates cross 
helicity conservation. A similar analysis [9j with hydro- 
magnetic load and metage yields conservation of helicity. 
The equation for an ideal non-dissipative MHD flow is: 

— + (v.V)v = ^ - VV> (39) 

at p pop 

where ip is the external potential and /io is the perme- 
ability of the medium. If s is the specific entropy and 
e(p, s) is the specific internal energy, one can write the 
following Lagrangian corresponding to the equation (|39|) : 



L = 



\v 2 -e{p,s) -ip 



B 2 
2p p 



pd 3 x (40) 
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The pressure P is related to e(p, s) via the thermody- 
namical relation: 



' -if 



The action is defined as: 



(41) 



They respectively suggest following variations inducing 
changes in p and B : 



Ap = -p(V.AX) (47) 
All, = BjdjAXi - BidjAXj (48) 



A = I Ldt 



(42) 



Following the notations developed in the section (III), 
as x is the position vector of a particle labelled by A = 
(Ai, A2, A3), the trajectories are given by: 



x = X(X,t) 



where X is a vectorial function that induces a mapping 
that if assumed one-to-one, can have an inverse. We 
would distinguish between the variations A incurred on 
a given label (i.e., at a point in the label space) and the 
variations 5 at a fixed point in the real space. If infinites- 
imal changes AX are made in the trajectories, for any 
quantity F one has the following relation [10(: 



SF = AF - AA.VF 



(44) 



Now, we have the continuity equation and the induction 
equation as follows respectively: 



Dp 
Dt 
DB 
~Dt 



dp 
dt 
dB 
~dt 



f (u.V)p = -p(V.v) (45) 
+ (v.V)B = (S.V)v - B(V.v) (46) 



Similarly, as is generally assumed if the specific entropy 
s of a fluid particle is constant, then 



As = 



(49) 



(43) 

Also, the corresponding variation of v is: 



Avi = —AXi 
Dt 



(50) 



Equipped with these relations, we proceed to apply vari- 
ation on the action (I42[) rewritten explicitly below 



A = / dt d x 
it-, Jv 



B 2 

~v z - e{p, s) - tp } p - — 



(51) 



We shall consider that suitably differentiable variations 
in the particle trajectories vanish for t outside [ti, £2] and 
AX.n = on the boundary of V. Therefore, we have 
starting from the expression (|5Tj) : 



5 A = 
5 A = 

SA = 



5A 



dt I d 3 x 



v ,U r-r~ r r / ! B - SB 
— — e — tp > op + p{v.ov - de — dip} 

2 J Mo 



(52) 



dt / d 3 x 
ti Jv 



^> j (Ap - AX. Vp) + p jtf.Atf - v(AX.V)v - Ae + AX. We - Atp + AX.Vip} 



5A 



in I. 



B.AB - B 



(AX.V)sj 



rff / d 3 x 



1 5£ r 

y - ^-(pe) -^j Ap + pv.Sv - p-^As - pVip.AX - AX.vj^y- - pe - pip 

If / -B 2 

^ B l {B J d ] )AX l - B l (B l d ] AX J ) - AX.V — 

Mo L V 2 



1 



V <i ^- - " W-} " " PVV - V - pe - pip )) .AX 



djiBiBj) + d t B 2 - \d t B z \ AXi 
Po I ^ 



d 3 xpv.AX 



d 3 Xv.AX 



t2 P 

dt I d 3 x 
tl Jt x Jv 



^(pv) + VP + pViP- ( -^}l + v(f- 
Dt po \ 2 A*o 



.AX 



(53) 
(54) 

(55) 
(56) 
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We pause here a bit and explain the steps involved in the 
above calculation. To reach from ([52"]) to ([53"]) , we have 
used ([34]): and ([53]) in turn used J48]) to yield (J54J). In 
arriving at ([55]), relations (gTJ), (gU) and ([10]) have been 
taken into account followed by integration by parts of the 
appropriate terms. The final step, relation ([56]) follows 
from (|55p by simple rearrangement and by using the rela- 
tion ([4"T]) . At this point, we argue that the variation is a 
symmetry displacement (i.e., a displacement that makes 
SA zero) and allow the equations ([3"9"]) and ([4"5]) to hold 
during the process of variation. Hence, the equation (l5(: 
yields: 



and therefore putting ([62]) in (|61[) . we arrive at the result: 



(v.AX)d 3 X — constant 



(57) 



Simply because one is always free to choose a new set 
of variables (A^, A' 2 , A3) such that 



d(A') _ d(X[,X' 2 ,X> 3 ) 
0(A) <9(Ai,A 2 ,A 3 ) 



(58) 



the A chosen to satisfy equation ([15]) is not unique. 
Hence, if the new set of variables does not modify the 
domain of integration, the value of the Lagrangian in the 
definition (|40|) remains the same. It is nothing but the 
fascinating concept of relabeling symmetry Q. 
As for the magnetic field lines lying entirely within the 
hydromagnetic fluid are always closed, we can translate 
MHD mctage /1 without affecting the Lagrangian; this 
means that the appropriate symmetry displacement is: 



AX 



dX 



A11 



(59) 



v.BdV — constant 



(63) 



In a rather more technical terms, one might choose to 
say that a one parameter translation group (a special 
subgroup of Arnold's symmetry group [l2j|) in the label 
space when represented by MHD load, MHD conjugate 
load and MHD metage has generated the conservation of 
cross helicity. 

VI. CONCLUSION 



In short, we have introduced the concept of magne- 
tohydrodynamic load in the case of MHD with untan- 
gled magnetic field lines and extended the definition for 
the non-trivial case of toroidal configuration of the mag- 
netic field lines in MHD. The elegant relation between the 
MHD loads and the magnetic helicity has been utilised 
to establish the relation ([37]) that very nicely connects 
the linking number of a pair of arbitrary magnetic field 
lines in the toroidal configuration discussed in the paper 
with the invariant loads in the magnetofluid having infi- 
nite conductivity. The MHD analogue of Ertel's theorem 
is rederived for the advected scalar - MHD load i.e., it 
has been shown that the potential magnetic field w.r.t. 
magnetohydrodynamic load is conserved. Moreover, the 
role of the MHD metage in generating the conservation 
of cross helicity in non- viscous non-resistive magnetofluid 
has been discussed. 



where, x — X(X,X,fi) is the transformation whose Jaco- 
bian is given by the inverse of the Jacobian given in the 
relation ([15]) . Thus, in the additional light of the fact 
that varying Lagrangian of the equation P0"]) by a con- 
stant A/x keeps L unaffected, the relation ([57]) becomes: 



OX q 

v.-^—Jd x = constant 
v d l l 

■* ( d *\ ,3 

v. p-r— \ d x = constant 
v \ 9 /i / 



But the relation (fl~4"]) suggests 11 
dX 



D 



(60) 
(61) 



(62) 



Acknowledgments 



One of the authors, SC, wants to thank Prof. J.K. 
Bhattacharjee for his constructive comments and discus- 
sions on an earlier draft of the paper. Also, the financial 
support from CSIR (India), in the form of fellowship to 
the same author, is gratefully acknowledged. Ayan's help 
in providing the authors with various relevant scholarly 
articles is heartily appreciated. 



[1] D. Lynden-Bell and J. Katz, Proc. R. Soc. Lond. A 378, Press, (1998) 

179(1981) [3] R. Salmon, Lectures on Geophysical Fluid Dynamics, Ox- 
[2] A.R. Choudhuri, The Physics of Fluids and Plasmas: An ford University Press, (2004) 

Introduction for Astrophysicists, Cambridge University [4] J. J. Moreau, C.r. hebd. Seanc. Acad. Sci. Paris 252, 



2810(1961) 

[5] H. K. Moffatt, J. Fluid Mech. 35, 117(1969) 
[6] A.F. Ranada, Eur. J. Phys. 13, 70(1992) 
[7] M. A. Berger and G. B. Field, J. Fluid Mech. 147, 
133(1984) 

[8] H. K. Moffatt and R. L. Ricca, Proc. R. Soc. Lond. A 



439, 411(1992) 
[9] A. Yahalom, J. Math. Phys. 36(3), 1324(1995) 
[10] F.P. Bretherton, J. Fluid Mech. 44, 19(1970) 
[11] E.A. Kuznetsov, JETP Lett. 76, 346(2002) 
[12] V.I. Arnold, J. Mec. 5, 19(1966) 



